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Compressible  flow  over  a  laser  turret  creates  a 
phase  distortion.  To  gain  insight  to  this  problem 
a  model  of  a  hemispherical  turret  is  developed.  The 
flow  field  is  obtained  using  a  second  order  solution 
based  on  Janzen-Rayleigh  expansion  for  the  compress- 
ible potential  flow  equation.  Contours  of  constant 
phase  distortion  were  calculated;  results  at  the 
critical  Mach  number  are  presented  and  discussed.  It 
was  found  that  the  distortion  at  0°  elevation  angle 
was  equivalent  to  a  positive  lens  whereas  at  90°  the 
distortion  is  similar  to  a  negative  lens.  At  4-5° 
the  main  effect  is  beam  tilting  into  the  wind.  For 
a  ratio  of  beam  radius  to  turret  radius  equal  to 
0.5,  the  focal  length  at  beam  elevation  angle  of  0°is 
16.5  km.  At  beam  elevation  of  54°  the  beam  tilt  was 
found  to  be  26.5  microradians .  The  results  of  this 
paper  are  compared  with  other  results  reported  in  the 
literature . 
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ABSTRACT 

Compressible  flow  over  a  laser  turret  creates  a  phase  distortion. 
To  gain  insight  to  this  problem  a  model  of  a  hemispherical  turret  is 
developed.  The  flow  field  is  obtained  using  a  second  order  solution 
based  on  Janzen-Rayleigh  expansion  for  the  compressible  potential  flow 
equation.  Contours  of  constant  phase  distortion  were  calculated;  results 
at  the  critical  Mach  number  are  presented  and  discussed.  It  was  found 
that  the  distortion  at  0°  elevation  angle  was  equivalent  to  a  positive 
lens  whereas  at  90°  the  distortion  is  similar  to  a  negative  lens.  At 
45°  the  main  effect  is  beam  tilting  into  the  wind.  For  a  ratio  of  beam 
radius  to  turret  radius  equal  to  0.5,  the  focal  length  at  beam  elevation 
angle  of  0°  is  16.5km.  At  beam  elevation  of  54°  the  beam  tilt  was  found 
to  be  26.5  microradians .  The  results  of  this  paper  are  compared  with 
other  results  reported  in  the  literature. 
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PHASE  DISTORTION  DUE  TO  AIRFLOW  OVER 
A  HEMISPHERICAL  LASER  TURRET 

by 

A.  E.  Fuhs*  and  S.  E.  Fuhs** 

I .  INTRODUCTION 

The  index  of  refraction  of  air  is  related  to  mass  density.  Starting 
at  modest  subsonic  Mach  numbers ,  e.g.  0.3,  and  extending  to  transonic 
flight,  flow  over  a  blunt  laser  turret  causes  significant  changes  in 
density.  The  change  in  density,  which  in  turn  changes  index  of  refraction, 
can  de focus  a  laser  beam. 

Variable  index  of  refraction  within  context  of  aeronautical  appli- 
cations is  not  a  new  problem;  however,  there  are  two  aspects  which  are 
new.   First,  the  laser  beam  is  large  relative  to  a  characteristic  turret 
dimension.   Second,  for  many  airborne  applications  the  laser  telescope 
requires  a  turret  which  protrudes  into  the  flow. 

Reconnaissance  aircraft  fly  photographic  missions.  The  compressible 
flow  over  the  aircraft  acts  like  a  distorted  lens  and  conceivably  could 
cause  distorted  photographic  images.  The  variable  index  of  refraction  is 
not  a  serious  problem  for  reconnaissance  aircraft  for  two  reasons.  The 
aperture  of  the  camera  is  small  compared  to  the  characteristic  scale  of 
the  distortions;  in  the  case  of  a  laser  telescope  this  in  not  true.  Fur- 
thermore, the  camera  on  a  reconnaissance  aircraft  is  not  mounted  in  a 
turret.  To  aim  the  camera  the  pilot  aims  the  aircraft.  Camera  windows 
are  flush  with  the  aircraft  surface. 

Cruise  missiles  with  intercontinental  ranges  may  use  celestial  navi- 
gation which  requires  precise  measurement  of  angular  location  of  stars. 


*   Distinguished  Professor  of  Mechanical  Engineering 
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A  distorted  lens  created  by  the  ambient  flow  field  causes  an  error  in 
measurement.  Knowledge  of  the  external  aerodynamics  is  important. 

A  model  which  is  useful  for  understanding  the  density  field  in 
the  vicinity  of  a  turret  is  the  flow  around  a  sphere  or  an  ellipsoid. 
The  plane  containing  an  incoming  fluid  velocity  vector  and  a  diameter 
of  the  sphere  does  not  have  flow  components  normal  to  the  plane.  This 
is  a  plane  of  symmetry  and  conceptually  can  be  considered  a  wall.  The 
flow  over  the  hemisphere  above  the  plane  can  be  calculated  readily. 

For  incompressible  flow,  the  solution  for  the  flow  over  a  sphere 

is  well  known.  See  for  example  pages  464-465  of  Milne-Thompson, 

(2)  .  (3) 

pages  92-93  of  Lamb    or  pages  339-342  of  Karamcheti.     Techniques 

are  available  to  introduce  the  influence  of  compressibility  or,  stated 

another  way,  the  influence  of  Mach  number.  One  method  is  the  Janzen- 

(4) 

Rayleigh  technique  discussed  briefly  by  Van  Dyke.     Another  brief 

(5) 

description  is  given  m  Chapter  10  of  Liepmann  and  Puckett.     The 

(fa) 

Janzen-Rayleigh  method  is  discussed  also  on  pages  328-334  of  Oswatitsch. 

(7) 

Lord  Rayleigh    provides  the  second  order  solution  for  a  sphere  in  a 

paper  published  in  1916. 


II.   DISCUSSION  OF  ANALYTICAL  MODELS 

Optical  distortion  due  to  the  external  flow  field  can  be  divided 
into  two  categories.  Viscous  flow  phenomena  fall  into  one  category  and 
include  shear  layers,  laminar  and  turbulent  boundary  layers,  and  shedding 
of  discrete  vortices.  The  other  category  involves  the  inviscid  flow 
field.  This  paper  treats  the  latter  category. 

Laser  turrets  can  be  classified  as  blunt  or  small-perturbation 
turrets.  If  the  turret  has  a  surface  with  a  normal  vector,  n,  aligned 
parallel,  or  nearly  so, to  the  freestream  velocity  vector,  V,  the  turret 
is  "blunt".   If  the  angle  between  the  normal  vectors  of  the  turret  sur- 
face and  the  freestream  velocity  vector  are  everywhere  large,  e.g.  60° 
to  90°,  then  the  turret  is  "small-perturbation". 

The  turret  being  considered  here  is  a  hemisphere,  which  is  obviously 
a  blunt  turret.  At  the  front  stagnation  point  of  a  hemispherical  turret, 
V  •  n  =  -  1.0. 

For  a  blunt  turret,  compressibility  effects  set  in  at  a  smaller 
freestream  Mach  number,  M^.  A  useful  concept  is  the  critical  Mach  number, 
M^".  At  the  critical  Mach  number,  somewhere  on  the  body  a  local  Mach 
number  is  just  sonic.  When  M^  >  M  *,  the  flow  becomes  inherently  nonlinear 
and  shock  waves  appear. 

For  a  circular  cylinder  with  its  axis  normal  to  the  flow,  the  critical 
Mach  number  is  0.3985.  For  a  sphere,  the  critical  Mach  number  is  0.5868. 
A  hemisphere-cylinder,  which  is  oriented  with  the  cylinder  axis  parallel 

to  the  freestream  flow  direction,  has  a  critical  Mach  number  less  than 

(8) 

0.7.  Hsieh    gives  experimentally  determined  pressure  coefficients  for 

a  hemisphere-cylinder  over  the  Mach  number  range  0.7  to  1.0.   In  refer- 
ence 9  Hsieh  extends  the  Mach  number  range  from  1.05  to  1.82.  Using  C 


from  reference  8,  the  maximum  local  Mach  number  is  1.25;  the  maximum 
local  Mach  number  occurs  at  an  angle  of  approximately  60°  on  the  sphere. 
The  angle  is  measured  from  the  stagnation  point  to  the  location  of 


minimum  C  . 
P 


For  small  perturbation  turrets  the  critical  Mach  number  is  larger. 
Calculations  were  made  of  a  turret  on  a  cylindrical  fuselage.      The 
turret  shape  was  determined  by 


R(x,6)  =  RQ  +  |  (1  +  cos  jt)    (1  +  cos  fe) 


R(x,6)  =  R  for  I  <      |x|   <  L  and  n/f  <  |e|<  tt 


(1) 


where  Rn  is  fuselage  radius;  e,  height  of  turret;  t,   length  of  turret; 
1/f,  the  fraction  of  the  2tt  circumference  occupied  by  the  turret;  and,  L, 
the  separation  between  periodic  turrets.  Table  I  summarizes  the  critical 
Mach  number. 


TABLE  I. 

Critical  Mach  Numbers  for  Small-Perturbation  Laser  Turrets 

l/R0 

t  An        n  is:                    /n         \t   •••       maximum  slope 
L/Rn        1/f       e/Rn        M  *           .     .  K 
0                  0         °°        on  turret 

1.005       5.0        .333 


0.05 

.88 

4.5° 

.10 

.82 

8.9° 

.15 

.76 

13.2° 

.20 

.72 

17.4° 

.25 

.68 

21.3° 

.30 

.65 

25.1° 

.35 

.62 

28.7° 

.40 

.59 

32.0° 

.45 

— 

35.1° 

.50 

— 

38.0° 

.55 

— 

40.7° 

There  are  different  analytical  models  that  can  be  developed  for 
laser  turret  geometries.   Flow  over  blunt  turrets  can  be  determined  using 
expansion  techniques  such  as  the  Janzen-Rayleigh  method  mentioned  ear- 
lier. As  more  terms  in  the  expansion  are  evaluated,  the  more  accurate 
are  the  flow  field  calculations  as  Mach  number  increases.  Upper  bound 
is  the  critical  Mach  number.   For  M  >  M  *,  the  calculations  will  be 

CO        CO   ' 

grossly  in  error. 

There  are  numerical  techniques  which  directly  integrate  appropriate 
differential  equations.  One  method,  which  is  applicable  to  axisymmetric 
flow,  is  the  technique  of  South  and  Jameson. 

For  small-perturbation  laser  turrets  the  equations  can  be  linear- 
ized. Compressibility  effects  in  subsonic  or  supersonic  flow  are  incor- 
porated. Transonic  flow  is  excluded  since  the  transonic  equations  are 
nonlinear  even  for  small  perturbation. 

Figure  1  summarizes  the  regions  in  a  body  slope  vs.  Mach  number 
map.  The  flow  within  the  region  ABCD  can  be  described  with  sufficient 
accuracy  using  the  linearized,  small-perturbation  equations.  The  line 
CD  represents  the  critical  Mach  number.  The  line  BC  is  shown  at  a 
slope  of  30° .  The  assumptions  involved  in  the  linearization  of  the  po- 
tential flow  equation  become  less  and  less  valid  as  the  maximum  slope 
increases.  The  region  CDJ,  which  resembles  an  inverted  triangle,  requires 
solution  of  the  nonlinear  small-perturbation,  transonic,  potential  equa- 
tion. 

The  region  EFGH  defines  the  area  where  Janzen-Rayleigh  technique  is 
most  useful.  The  line  GH, which  defines  the  critical  Mach  number,  is 
shown  dashed;  the  critical  Mach  number  is  a  function  of  turret  shape. 
Hence  GH  is  intended  to  qualitatively  suggest  the  upper  bound  for  Mach 
number. 

5 


3  4 


ill 
a. 
O 


90 
80 


</>    60 


>- 
o 
o 
co 


3 
2 

X 

< 


40 


20 


T\ 


i 


i    i 


BLUNT 
TURRETS 


i 


B  I  2  34 

SMALL-PERTURBATION 
TURRETS 

A I i L 


0.2  0.4  0.6  0,8  1.0 

FREESTREAM  MACH  NUMBER  ,    M^ 


Figure  1.      Freestream  Mach  Number  Regions  for  Application  of 
Various  Analytical  Models. 
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There  are  vertical  lines  1,  2,  3,  and  4  within  region  EFGH.  These 
lines  represent  the  upper  bound  of  applicability  of  the  first,  second, 
third,  and  fourth  order  solutions  when  Janzen-Rayleigh  technique  is 
used.  The  first  order  solution  may  be  applied  within  region  F11E;  the 
second  order  solution,  within  region  F22E;  etc. 

A  hemispherical  turret  considered  in  this  paper  is  an  example  of 
a  blunt  turret  falling  along  the  line  FG. 


III.  JANZEN-RAYLEIGH  EXPANSION  TECHNIQUE 
In  spherical  coordinates,  the  equations  of  motion  are  as  follows 

(2) 


2  2 

3u  JL  v  3u   v       1  3d     a.   9p 


3r   r  30   r      p  3r     p   3r 


2 
^+v_3j/  +  uv__ll3p_a_3£  (    . 

3r   r  36   r      p  r  39     pr  39 


1  f  3p  ,  v  3p"|     13.2s     1     3  .        .   Q.   n   ,„» 
p-  [d  3?  +  v     Je\     +     ^     3F  (r  U)  +  FiETe  3?  (v  Sin  9)  =  °   (4) 

Equations  (2)  and  (3)  are  Euler's  equations.  Equation  (4)  is  the  con- 
tinuity equation.  For  flow  over  the  sphere,  there  is  no  dependence  on 
the  azimuthal  angle  in  polar  coordinates.  Consequently,  equations  (2) 

and  (4)  are  two-dimensional  equations  depending  only  on  r  and  9. 

2  2 

Multiplying  equation  (2)  by  +  u/a  and  equation  (3)  by  -  v/a  , 

adding  and  combining  with  equation  (4)  yields 

13,2.     1      9  .         .   Qx    1  )    2   3u  L  uv  3v  ,     3v  x  v  3v   ,n 

— 2  —  (r  u)  +  . — -  —  (v  sin  9)  =-2<u   —  +  —  -tt^+uv— -  +  —  —  >  (5) 

rz   3r        r  sin  9  39  az  }    3r   r  39      3r   r  39  j 

The  speed  of  sound,  a,  is  the  local  value  and  changes  from  point  to  point 
in  the  flow.  It  is  necessary  to  account  for  variations  in  a.  Using  the 
energy  equation  along  a  streamtube,  one  can  show  that 

4  =  i  +  *£  m2  (i  -  ni^-Z  )  (6) 

a2        2    °°        U2 

Inviscid  flow  is  considered.  A  potential  can  be  introduced  for  the 

velocity  components  ,,  _3J>_  (7\ 

3r 


where  U  is  the  freestream  velocity  at  upstream  infinity.   Introducing 

equations  (7)  and  (8)  into  equations  (5)  and  (6)  yields  equations  in 

terms  of  a  potential. 

The  Rayleigh-Janzen  expansion  considers  the  potential  function  to 

be  given  as 

<f>  =  4>n  +     M2     <K    +  M4  cf>0  +  (9) 

T        U  ooTJ_        °°Z      ... 

The  solution  for  ^  is  considered  the  first  order  solution.  The  solution 
involving  both  <J>n  and  <j>..  is  called  the  second  order  solution.  In  this 
paper  only  the  second  order  solution  is  considered.  For  ease  of  writing, 
the  following  definition  is  introduced: 

\l>   =  M2  <j>_  (10) 

00     J_ 

To  digress  briefly  from  the  development  of  the  Janzen-Rayleigh 

expansion,  the  solution  for  the  flow  over  a  cylinder  normal  to  the  fluid 

(<4)  . 
velocity  will  be  discussed.  Van  Dyke    indicates  the  solution  for  the 

cylinder  has  been  performed  to  M^  .  For  y   =  7/5,  the  result  for  maximum 

surface  speed  is 

^lax  ?  4 

k  =  2.00000  +  1.16667M  +  2.57833M  +  


U 

+  25.59041M  +....+  7450M16 


(11) 


The  critical  Mach  number  for  a  cylinder  is  Mto  =  0.39852  +  0.00020. 

In  terms  of  the  potential  function,  the  speed  of  sound  becomes 

a2        2 

2g-  =  1  +  K   FC6  ,  *  ,  r,  Y)  (12a) 

a 

00 

where  F  is  given  by 


F  = 


Y-l 

2 


[-  a  4)] 


(12b) 


With  substitution  of  equations  ( 9 )  and  ( 10 ) ,  the  potential  function 
becomes 


V2(f)n  +  V2^  =  M2 
0  °° 


9  2  1     o  1 

r  rr   rz  r  8  Or   r^  6  r   r 


(13) 


where  all  terms  on  the  right-hand  side  are  determined  from  the  first 
order  potential  function  <)>„.  The  result  of  equating  coefficients  of  like 


powers  in  M^  yields  two  equations 

V2d>n  =  0 


and 


V2ii>  =  M2  V2d>     =  M2 

OO  -I  oo 


0 


o  2  1 

r  rr       r       r  0   0  r       r 


,2    d>    +  — u    d>2    d> 

6  ^r    r4  Te  y6f 


(14) 


(15) 


The  right-hand  side  of  equation  (15)  is  known  due  to  the  fact  that 
equation  (14)  can  be  solved  to  yield  d>n.  To  emphasize  this  fact,  the 
first  term  on  the  right-hand  side  is  written  explicitly  showing  <f>n, 


4>2<i> 
r  rr 


m 


924>( 
dr77 


Equation  (15)  is  Poisson's  equation. 


IV.   FIRST  ORDER  SOLUTION 


The  first  order  solution  is  obtained  from  Milne-Thompson,    Lamb, 


or  Karamcheti. 


(3) 


It  is 


f  AN      TT  /  Q   L  R3   COS   9   N 

>Q(r,  9)  =  U  (r  cos  6  +  — ^ ) 


(16) 


where  R  is  the  sphere  radius.  The  upstream  stagnation  point,  P,  is  at 
P(R,  180°);  and  the  downstream  stagnation  point  is  at  P(R,  0°).  Using 
equation  (16)  one  can  determine  the  derivatives  to  substitute  into  equa- 
tion (15).  These  are 


R3 

)  =  u  cos  e  (l  -  £*  ) 

r  rJ 


rr 


=  U  cos 


3R3 


r  U  sin  0 


r  U  cos  6 


1  + 


1  + 


1 

R3~ 

2 

r3 

1 

2 

R3" 

0r 


=  -  U  sin 


1  - 


R; 


(17) 

(18) 
(19) 

(20) 
(21) 


The  result  of  combining  equations  (15)  and  (17)  to  (21)  is 


V2,  =  £2LA 
1     r 


J3A(1-A)2  +  sin2  6 


4  -  ~   A  +  6X2 


-¥■*]) 


(22) 


,3,  3 


where  X  =   R  /r  has  been  introduced. 

Before  proceeding  with  second  order  solution,  a  comment  about 
boundary  conditions  is  appropriate.  The  component  of  velocity  normal  to 
the  surface  of  the  sphere  must  be  zero.  At  infinity,  the  velocity  should 
be  U.  Consequently, 


3<j> 


1 


3r 


R 


=  0  ; 


'  r  88 


=  0 


(23) 
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Since  <j>Q  yields  U  at  infinity,  the  velocity  perturbation  due  to  <£, 
should  be  zero  at  infinity. 
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V.  RAYLEIGH  SOLUTION  FOR  SECOND  ORDER  FUNCTION 

(7) 

In  July,  1916,  Lord  Rayleigh    reported  the  solution  for  the 

second  order  function. 

The  Legendre  functions 

P,  =  cos  e  (24) 

P3  =  |  cos3  e  -  |  cos  6  (25) 

are  introduced  into  equation  (15)  in  the  following  manner:  First, 
equations  (17)  through  (21)  are  substituted  into  equation  (15);  second, 
terms  are  grouped  so  as  to  form  the  functions  P-,  and  P^.  The  result  is 


^  -  U 


(     36  R6       9  R9      \  /     R3       24  R6        3     R9    \        " 

\~  5~  W  +  2  ~W)  Pl       [b  FT  "  5~  rT  +  2"    rWJ    P3 


(26) 

IX  II"*  0        1"  Z        I'l  u     I  0  1 

Equation  (26)   is  in  the  form 


V2cj>,    =  Z     £     r      PC  (27) 

1  n  nm 

n    m 

where  C.  are  zero  except  for  m  =  7  and  10;  all  C9  are  zero;  and  C~ 

are  zero  except  for  m  =  4,  7,  and  10.  Non-zero  C   are  the  appropriate 

constants  appearing  in  equation  (26).  The  Legendre  differential  equation 

is 


d  dP 


sin  6  de 


[Of    -i 
sin  6  dTj  +  n  (n+1)  Pn  =  ° 


(28) 


Define  the  function  d>   which  satisfies  the  equation 

nm  M 

324  „     3d)  ,  T  94>     1 

-nu  r       _     2x       _         n™  j    2       Tnm  ,       1  3      I    .      .     Tnm  /onN 

r   Vnm  "  v  V  "  3FT-  +  ?   W~  +  r^ilFT  ae  [sin  9  ae — J         (29) 

Assume 

=  r2"m  P    A  (30) 


nm  n    nm 
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Consequently, 


3c}> 


nm  _  ,  0  x  1-m  _  . 
r— -  =  (2-m)  r  PA 
3r  n  nm 


(31) 


and 


324> 


nm  _  ,0  ,  ,  .  -m  „  . 
r— 5—  =  (2-m)  (1-m)  r  P  A 
3r2  n     nm 


(32) 


Combining  equations   (29)  and  (32)   gives 


r_n¥  C   =  (2-m)  (1-m)  r  m  P  A 


n  nm 


n  nm 


(33) 


+  2  (2-m)  r  mP  A   +  2 


-m+2 


[ 


3P 


sin  6 


n 


A 


nm 


n  nm   r2  sin  6  99        39 
Substitution  of  equation  (28)  into  equation  (33)  yields 


•»[■ 


rmPC   =  A   rmP|  (2-m)  (1-m)  ■■  2  (2-m)  -  n  (n+1) 

n  nm    nm 


■] 


(34) 


As  a  result 


C 


A 


nm 


nm 


(2-m)  (3-m)  -  n  (n+1) 


(35) 


Using  equations  (26),  (30),  and  (35),  that  part  of  the  second  order 
function  identified  as  the  particular  solution  is  obtained 


\   =  U 


h 


R6  +  R9 


]•  '{• 


3R3    3R6    3R9 
10r2   I0r5   I76F" 


5r5   ^r8 
The  homogeneous  solution  is 


(t>1  =  Ar  P     +  BP1/r2  +  Cr  P3  +  D  P3/r1+ 


(36) 


(37) 


At  infinity  the  potential  must  be  finite  and  the  velocity  must  be  U; 
the  function  <f>n  of  equation  (16)  yields  U.  Consequently  both  A  and  C 
equal  zero.  At  the  surface  of  the  sphere 
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3<|> 


1 


i   -0 
(r=R 


and 


r  36 


=0 


(38) 


(9=0 


From  the  requirements  of  equation  (38)  combined  with  equations  (36) 
and  (37) 


R3 

B  =  §-  U  M2 
3     °° 


and 


D  = 


27R5U  m2 
55    °° 


The  final  expression  for  the  potential  function  is 


'0+M-*l=fU 


r   +  *-z 


R3 

2r' 


Ll 


+  UM 


R3    R6  +  R9 


3r2   5r5   24r8 


3R; 


27R5   3R6    3R9 


10r2   55?^   10r5   176r8 


(39) 


(40) 


(41) 


The  first  term  in  brackets  is  <j)n;  and  the  second,  long  term  is  the  second 


0 


order  contribution. 
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VI.  DEVELOPMENT  OF  PHASE  DISTORTION 


Use  of  Potential  Function  to  Obtain  Index  of  Refraction 
The  index  of  refraction  and  density  are  related  by 


n  =  1  +  k'  ^—  =  1  +  K  £_ 


P«  p 


SL 


-4 


(42) 


where  «'  is  a  constant  with  value  of  about  2  x  10  .  Note 

k  =  k'p^/pot  •  The  flow  is  assumed  to  be  isentropic.  Consequently 


1  +  l£  M2 

2   °° 


1  +^M2 


Y-l 


(43) 


The  local  Mach  number  is  given  by 


M2  =  i 


lary   yr  ae y 


(44) 


In  equation  (42)  p  is  the  local  density,  p  is  the  freestream  density, 
and  pCT  is  the  sea  level  density  on  a  standard  day.  By  combining  equa- 
tions  (6),  (12a),  (41),  (42),  (43),  and  (44),  an  expression  is  obtained 
for  the  index  of  refraction. 

The  turret  has  a  characteristic  length;  for  the  sphere  this  is 
the  radius  R.  The  laser  beam  has  a  characteristic  size  which  is  the 
diameter  D  or  radius  R' .   The  wavelength  of  the  laser  radiation,  A, 
is  also  a  characteristic  length  as  will  be  demonstrated  in  subsequent 
discussion.  Two  nondimensional  lengths  can  be  formed.  Since  R  occurs 
prominently  in  the  equations  for  the  potential,  R  is  chosen  as  the 
reference  length.  With  that  choice,  A/R  and  D/R  are  two  of  the  para- 
meters specifying  the  problem. 
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Phase  distortion  is  a  function  of  the  direction  that  the  beam  is 
pointed.  Define  £  and  a  as  elevation  and  azimuth  angles.  For  the 
hemispherical  turret,  symmetry  exists  about  the  diameter  of  the  sphere 
which  is  parallel  to  the  freestream  velocity  vector.   In  general,  two 
angles,  e  and  a,  are  needed  to  define  the  beam  direction •  however, 
for  the  hemispherical  turret  only  one  angle  is  needed.  Looking  ahead, 
the  angle  6'  occurring  in  equations  (45)  and  (46)  will  be  used  as  the 
beam  angle;  see  also  Figure  2. 

Geometrical  Relations  for  Calculation  of  Phase  Distortion 

Figure  2  shows  the  geometry  of  the  problem.  The  radius  of  the 

sphere  is  R.  The  beam  is  elevated  at  an  angle  6'.  The  z  axis  is  the 

same  as  the  polar  axis  for  polar  coordinates. 

To  locate  points  within  the  beam,  three  variables  are  used;  s  is 

distance  along  the  beam,  as  shown  in  Figure  3;  a  is  an  angle  measured 

clockwise  from  the  windward  edge  of  beam;  and  R'  is  the  radial  distance 

from  the  beam  centerline  or  axis. 

The  flow  properties,  i.e.,  density  and  local  Mach  number,  are 

functions  of  r  and  6.  These  are  polar  coordinates.  Consequently,  it  is 

necessary  to  express  r  and  0  in  terms  of  6',  s,  a,  and  R.  The  equations 

are     2   r  l2     2     2 

r  =   (R  +  s)  sin  6'  -  R'  cos  a  cos  0'   +  R'   sin  a 


2  (45) 


and 


.   -1 
=  tan 


=  (R  +  s)  sin  0'  -  R'  cos  a  cos  0 'J   +  R' 
+  (R  +  s)  cos  0'  +  R'  cos  a  sin  6' 

VI  2     2     2 
(R  +  s)  sin  0'  -  R'  cos  a  cos  0'    +  R'   sin  a 


(46) 


(R  +  s)  cos  0T  +  R'  cos  a  sin  0' 


Equations  (45)  and  (46)  can  be  used  to  evaluate  equations  (41)  to  (44). 
The  index  of  refraction  is  given  by  equation  (42).  Define  N  as  the 
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BEAM 


{j>FLOW 


Figure  2.  Geometry  for  Phase  Distortion  Calculation. 


Figure  3.  Cross  Section  of  the  Turret  and  Beam  in  the  x-z  Plane. 
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number  of  waves  along  a  distance  s.  Hence 

dN  =  p-  (47) 

A 

where  A  is  the  local  wavelength.  The  local  wavelength  is  a  function  of 
0' ,  a,  R' ,  and  s.   It  is 

A  (1  +  k) 

A  "  f  "  Hf  "  f(l  +  KP/p  )  "  (1  +  Kp/p  )  ^^ 

oo  oo 

As  a  result,  the  number  of  waves  is 


/ 


:  x_xi\  K)  (1  +  K  r> 


-s 

A  reference  is  needed.  The  reference  which  is  convenient  is  the  number 

of  waves  along  the  beam  axis;  designate  that  quantity  as  1\L.  Using  this 

reference 

s  0  s 

ft.   v  /  /  (l+<  —  )ds  -  /  (l+K^)ds  -  /  (1+k— )ds 
0   A  C1+k)  1  /       p        /      p       fp 


Note  the  density  on  the  beam  centerline  for  -s  '<_  s  <_  0  is  denoted  by  p  , 
a  reference  density.  For  0  <_  s  <  °°,  the  density  on  the  centerline  is 
Pn  which  is  a  function  of  s.  In  equation  (50)  the  value  of  unity  in 
each  integrand  can  be  cancelled  giving 


N"No  =  rrfer- 


The  second  integral  in  equation  (51)  has  been  termed  the  GAP  INTEGRAL 
since  it  enters  the  equation  as  a  result  of  the  gap  between  s  =  -s'  and 
s  =  0 .  See  Figure  3 .  In  the  calculations  which  were  conducted ,  the  gap 
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integral  was  evaluated  by  assuming  p  =  p  ;  this  is  an  arbitrary  assump- 
tion. Note  that  in  equation  (51)  wavelength  is  nondimensionalized  by  R 
to  form  X^/R.  The  distance  along  a  ray  is  also  nondimensional ,  i.e., 
ds/R. 

A  computer  program  has  been  developed  for  the  HP9830  which  calculates 
N  -  1\L  as  a  function  of  sphere  radius,  R;  elevation  angle,  8';  radius 
within  the  beam,  R' ;  angle  within  the  beam,  a;  freestream  Mach  number, 
M;  wavelength,  X^;  constant  in  density  equation  (42),  k';  ratio  of  heat 
capacities  for  air,  yj  and  freestream  speed  of  sound,  a^.  Note  that 
it  is  not  necessary  to  specify  p^  since  all  calculations  depend  on  den- 
sity ratio;  see  equations  (48)  and  (51).  Values  used  in  the  calculations 

are  as  follows: 

R  =  0.9144  meter 

— fi 
X     =3.8  micron  =   3.8  x  10       meter 

oo 

k'    =   2  x  10"4 

a  =  342  meter/ sec 

Y  =  1.4 
M  ,  R' ,  a,  and  6'  have  been  varied.  Calculations  have  been  performed 
to  the  extent  necessary  to  plot  phase  distortion  maps. 


19 


VII.  RESULTS  OF  CALCULATIONS 

Graphical  Presentation  of  Results 

Plots  of  isocontours  of  phase  shift  were  made  for  steps  of  18° 
starting  with  6'  =  0  to  9'  =  90°.  The  plots  are  shown  in  Figures  4 
to  9.  To  avoid  awkward  decimal  values,  units  of  the  plots  are  deci- 
wavelength;  hence  250  is  2.5  wavelengths,  and  25  is  a  quarter  of  a 
wavelength.  Recall  that  the  beam  center  line  is  used  as  reference  for 


Figure  4.  Contours  of  Constant 
Phase  Shift  for 
6'  =  0°. 


Figure  5.  Contours  of  Constant 
Phase  Shift  for 
9  =  18°. 


the  phase  shift.  Reference  to  equation  (51)  shows  that  a  positive 
phase  shift  is  the  case  of  more  waves  along  the  ray  in  question  than 
along  the  ray  coincident  with  the  beam  centerline. 

For  the  series  of  plots  in  Figures  4  to  9  the  critical  Mach 
number  was  selected.  For  a  sphere  this  is  M*  =  0.587.  The  gap  inte- 
gral is  not  included  in  the  phase  shift.  The  outer  edge  of  each  plot 
is  for  D/2R  =0.5.  Recall  that  D  =  beam  diameter,  and  R  =  turret  radius. 
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Figure  6.  Contours  of  Constant 
Phase  Shift  for 
9'  =  36°. 


Figure  7.  Contours  of  Constant 
Phase  Shift  for 
e!  =  54°. 


220' 


Figure  8.  Contours  of  Constant 
Phase  Shift  for 
0*  =  72°. 


Figure  9.  Contours  of  Constant 
Phase  Shift  for 
6'  =  90°. 
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To  provide  insight  to  the  magnitude  of  gap  integral,  Table  II 
gives  the  maximum  value  for  each  of  the  plots. 

TABLE  II.  MAXIMUM  VALUE  OF  GAP  INTEGRAL 


beam  elevation  angle 


maximum  value  of  gap  integral 
in  deciwave  lengths  and  location 
within  beam 
(gap  integral) a 


0° 
18° 
36° 
54° 
72° 
90° 


+  36.8 
+  65.4 
+  67.8 

-  83.6 

-  86.3 

-  86.6 


0° 

0° 

180° 

120° 

90° 


M  =0.586  for  these  calculations.  Also  D/2R  =  0.5. 

oo 

"  not  a  function  of  a. 

Evaluation  of  Phase  Distortion  as  a  Function  of  Distance  Along  the  Beam 
An  important  question  is  the  rate  of  decay  of  the  integrand  in 

equation  (51).  Calculations  were  made  of 

3(N  -  NQ) 


a(s/R) 


which  is  the  phase  distortion  per  unit  distance  along  the  beam.  The 
phase  distortion  per  unit  distance  is  a  function  of  distance  along  the 
beam.  Results  of  the  calculation  are  shown  in  Figure  10.  The  phase 
distortion  is  the  area  enclosed  by  one  of  the  curves.  Most  of  the  phase 
distortion  occurs  within  a  length  equal  to  one  turret  radius. 
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Figure  10 . 


0.25        0.50        0.75  1.0  1.25 

RATIO  OF  DISPLACEMENT  AL0N6   BEAM 
TO  TURRET    RADIUS,      S/R 

Phase  Distortion  Per  Unit  (s/R)  as  a  Function 
of  Distance  Along  the  Beam  for  Rays  at 
a  =  0°,  60°,  120°,  and  180°  Within  the  Beam. 
Elevation  Angle  was  9'  =18°. 


Dependence  of  Phase  Distortion  on  Flight  Mach  Number 

Another  parameter  of  considerable  interest  is  the  dependence  of 
phase  distortion  on  Mach  number.  Qualitatively,  isodensity  contours  at 
M^  <  M*  have  an  appearance  the  same  as  for  Moo=  M*.  Quantitatively,  there 
are  changes  in  both  the  shape  of  an  isocontour  of  constant  phase  shift 
and  the  magnitude  of  the  phase  shift.  Table  III  compares  magnitudes 
of  distinctive  features  of  the  maps  for  two  Mach  numbers.  Table  III 
should  be  studied  simultaneously  with  the  appropriate  figure  for  0 ' . 

Phase  distortion  curves  based  solely  on  a  first  order  calculation, 
i.e.,  using  <f>n  in  equation  (9)  only,  should  have  the  same  shape  for  all 
M  .  A  calculation  of  phase  distortion  based  on  a  second  order  solution, 
i.e.,  using  both  <j>n  and  <i>,  in  equation  (9),  should  be  a  function  of  M^. 
Figure  11  shows  the  change  in  shape  of  a  phase  distortion  curve  with 
Mach  number.  The  curves  are  drawn  for  6'  =  18°;  M  ,  either  0.300  or 
0.587;  R'/R  from  0  to  0.5  and  for  a  =  0°  and  a  =  180°.  The  phase  dis- 
tortion has  been  normalized  with  the  value  at  R'/R  =  0.5  and  a  =  180°. 
Differences  in  the  curves  are  due  to  Mach  number  effects. 
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TABLE  III.  COMPARISON  OF  MAGNITUDE  OF  PHASE  SHIFT  FOR  TWO 
DIFFERENT  MACH  NUMBERS 

feature  of  isocontour  map     value  at  M  =  0.300   M  =  0.587 

*  on  oo 


0° 
18° 

36° 

54° 

72° 

90° 


maximum  <J>  at  beam  edge 

maximum  c|>  at  positive  peak 

minimum  <$> 

maximum  <J> 

minimum  $ 

maximum  <J> 

minimum  $ 

maximum  <J> 

minimum  <f>  at  negative  valley 

maximum  tf>  at  a  =  0° 


-41.5 

+15 

-69.6 

+44.2 

-80.3 

+85.2 

-54.4 

+  83.7 

-14.1 

+49.1 


-160 

+  60 

-290 

+240 

-345 

+370 

-268 

+386 

-   74.1 

+  223 


M.«  0.300 


NONDIMENSIONAL  RADIUS 


Figure  11.      Phase  Distortion  for  Two  Mach  Numbers, 
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Interpretation  of  Results 

Phase  shift  for  9 '  =  0  as  shown  in  Figure  4  consists  of  a  series 
of  concentric  circles.  Since  the  phase  shift  is  negative,  the  beam 
is  being  focused.  Concentric  curves  with  negative  phase  distortion 
gives  a  positive  focus. 

The  contours  of  constant  phase  shift  for  9 '  =  54°  are  almost 
straight  lines.  The  algebraic  signs  of  phase  shift  contours  indicate 
the  beam  is  being  tilted  in  a  direction  opposite  to  the  relative  wind. 
The  angular  tilt  is  given  roughly  by 


(tilt  angle)  =  KN-Nq)^  -  (N-N0)180O"| 


A 

(52) 


2R1" 

Inserting  the  appropriate  values  (values  for  N-Nn  were  obtained  from 
Table  III)  gives 

...nj.    ,  x    (370  +  268)  (3.8  x  10_6m)   oc  c  •      ,. 
(tilt  angle)  =  (2)  (.4572m)  (100) =  26.5  microradians 

Tilt  of  a  beam  does  not  cause  a  serious  problem. 


Cook 


These  results  for  tilt  angle  can  be  compared  with  similar  work  by 
(12) 

(12) 
Captain  Cook    used  the  following  values  m  a  sample  calculation 


(notation  has  been  adjusted  to  match  this  report) : 

R'  =  0.25m  X     =10.6  micron 

00 

R  =  0.71m  k'/p^  =  0.22  cm3/gm 

M  =  0.5  k'=2.7x  10~4 

00 

For  9'  =  54°,  a  beam  tilt  of  63  microradians  was  calculated.  When  adjust- 
ments are  made  for  different  parameters  used  to  calculate  the  two  results, 
i.e.,  63  and  26.5  microradians,  the  answers  agree  quite  favorably. 

It  is  possible  to  determine  an  approximate  focal  length,  F  .  The 
subscript,  a,  indicates  the  focal  length  in  a  plane  defined  by  a  of 
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Figure  2.  The  appropriate  equation  is 


TWJ2 


(N-NJ   -  2(N-Nn) 


0 


R2 


0' 


Ri 


(53) 


Equation  (53)  was  evaluated  for  a  =  0  and  for  the  following  values  for 

other  parameters 

-3      3 
p  =  1.23  x  10   gm/cm 


A  =  3.8  x  10   m 


AR'=  0.1524  m 


M  =  0.587 


R  =  0.9144  m 


(N-NQ)  at  R^  =0.799  wavelengths 

(N-NQ)  at  Rj  =0.214  wavelengths 

9f  =  0° 
A  focal  length  of  16.5  km  was  found. 

Before  proceeding,  some  additional  information  concerning  equation 
(53)  is  appropriate.  Equation  (53)  is  a  specialization  of  the  equation 


x 


3x2 


(54) 


where 


b 


ds 


(55) 


For  the  problem  at  hand  equations  (54)  and  (55)  yield 


i    <     (  rp2-po  .     .  r  prpo  .    r 


p0"p0 


ds 


(56) 


where  p„  =  p(x?),  p,  =  p(x-,),  etc.  It  is  obvious  that  the  third  integral 

(12) 
is  zero.  Equation  (56)  leads  to  equation  (53).  From  Cook's  "   report  for 

the  previously  cited  values  for  R,  R' ,  M^,  etc. ,  a  focal  length  of  about 

6  km  is  calculated.  Adjusting  the  value  for  values  of  R,  R' ,  M  ,  etc., 
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used  in  the  calculations  of  this  report,  good  qualitative  agreement 

is  obtained. 

(13) 

Wolters  and  Laffay     have  calculated  focal  length  for  an 

aperture  of  a  turret  on  a  F-15.  The  turret  geometry  was  a  cylinder 
with  axis  parallel  to  the  freestream;  radius  of  the  cylinder  was  0.5m. 
When  viewed  from  a  position  normal  to  the  axis,  a  cylinder  appears  as  a 
rectangle.  The  square  corners  were  rounded  with  a  radius  of  0.3  m. 
Without  rounding  of  the  corners,  the  cylinder  viewed  from  its  end  would 
appear  as  a  disc  with  radius  0.5  m.  With  the  rounding  of  the  corners, 
the  disc  is  reduced  to  0.2  m;  the  disc  is  normal  to  the  freestream 

velocity  vector. 

(13) 

Wolters  and  Laffay   '  determined  the  focal  length  of  the  lens 

created  by  density  variations  within  the  external  flow.  The  following 
parameters  were  used  in  their  calculations: 

R'  =  0.17  m  X  =  10.6  x  10~6  m 

00 

R  =  0.5  m  k'/p  =  0.223  x  10"3  m3/kg 

si  00  ° 

M  =  0.587  9'  =  0 

00 

altitude:   20,000  ft. 
For  these  conditions  a  focal  length  of  18.2  km  was  determined.  Since 
the  distortion  varies  linearly  with  the  ambient  density,  a  fact  which 
follows  from  equation  (42),  the  focal  length  can  be  adjusted  to  sea  level 
by  an  adjustment  to  density.  For  the  same  conditions  as  above  except  at 
sea  level,  the  focal  length  would  be  8.2  km.  A  value  8.2  km  is  consistent 
with  results  obtained  in  this  report,  16.5  km. 

Having  compared  in  a  rough  way  calculated  results  from  references 
(12)  and  (13)  with  calculated  results  from  this  paper,  the  relative 
advantages  and  disadvantages  of  each  model  will  now  be  compared.  Reference 
(12)  used  al/3x  and  3l/3y  to  predict  tilt;  92I/3x2  [see  equation  (54)]  and 
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32I/3y2  were  used  to  predict  focal  length.  These  were  evaluated  at  the 
beam  center.  An  analytical  solution  was  obtained.  Reference  (12)  also 
used  a  first  order  solution  for  the  flow  field,  which  is  equation  (16). 
This  paper  uses  a  second  order  solution  as  given  by  equation  (41).  Hence 
a  more  accurate  flow  field  is  used  in  this  paper.  Reference  (12)  yields 
a  solution  in  simple  form;  this  paper  requires  a  desk  top  computer,  eg. 
HP9830,  or  else  a  high  degree  of  patience  if  a  HP55  were  to  be  used. 

In  contrast  to  reference  (13)  which  used  a  numerical  solution  to 
the  flow  field,  the  second  order  analytical  solution  gives  the  flow  field 
to  any  degree  of  spatial  resolution  required.  Numerical  solutions  give 
results  in  discrete  steps.  Based  on  Figure  5  of  reference  (13)  the  ratio 
of  step  size  to  cylinder  radius  is  0.22.  Figure  10  of  this  paper  shows 
that  approximately  one-half  of  the  phase  distortion  is  generated  in  a 
step  size  0.22.  Of  course  a  smaller  step  size  can  be  used.  Numerical 
techniques  have  the  advantages  of  providing  solutions  for  complex  turret 
geometry.  This  paper  is  applicable  only  to  a  hemispherical  turret. 

Look  at  Figure  9  for  which  6'  =  90°.  The  contours  have  a  symmetry 
about  the  a  =  0°--180°  line  and  about  the  a  =  90°- -270°  line.  The  phase 
distortion  has  positive  values  which  give  a  negative  lens.  The  focal 
length  for  a  =  0°  is  considerably  less  than  for  a  =  90°.  At  6 '  =  90°, 
the  beam  is  defocused. 

It  is  apparent  from  Figures  4  to  9  that  the  phase  distortion  is  a 
strong  function  of  elevation  angle  0 ' .  Compensation  for  the  phase  dis- 
tortion is  possible,  at  least  partially.  However,  the  shape  of  the  pri- 
mary telescope  mirror  required  for  compensation  is  a  function  of  both 
elevation  angle  and  Mach  number.   Such  compensation  should  fall  within 
the  capability  of  adaptive  optics  techniques. 
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To  discuss  the  altitude  dependence,  recall  <  =  k'p^/Pot  •  Both 
<'  and  Pot  are  constants  for  a  given  wavelength  and  specified  atmospheric 
conditions.  Consequently,  k  varies  directly  with  p  .  Now  examine  equa- 
tion (51).   Since  <  <<  1,  the  term  1  +  <  is  essentially  unity;  it  follows 
that  N-Nn  varies  directly  as  k.  Phase  distortion  changes  as  p^.  Obviously 
as  altitude  increases,  phase  distortion  decreases. 
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VIII.  CONCLUDING  REMARKS 

Compressible  flow  over  a  sphere  has  been  obtained  forming  the  basis 
for  a  model  of  a  blunt  hemispherical  laser  turret.  Phase  distortion  has 
been  calculated  for  several  freestream  Mach  numbers.  Phase  distortion 
data  are  presented  for  M  equal  to  the  critical  Mach  number.  Calculations 
are  consistent  with  those  of  references  (12)  and  (13). 

The  phase  distortion  is  a  function  of  both  elevation  angle  and  Mach 
number.  Hence  optical  compensation  requires  a  variable  shape  primary 
mirror. 

The  model  of  a  hemispherical  laser  turret  provides  insight  to  the 
behavior  of  phase  distortion  as  a  function  of  Mach  number.  See  Figure 
11.  Also  refer  to  Table  III. 

Another  valuable  aspect  of  the  calculations  is  the  decay  of  den- 
sity perturbation  with  distance  along  the  beam.  See  Figure  10.  It  is 
apparent  that  the  distortion  is  fairly  localized.  Within  a  distance  of 
s/R  =  1.0,  over  99%  of  the  distortion  occurs. 

The  formulation  clearly  shows  the  altitude  dependence  of  phase 
distortion. 
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